The dissipative two-level system under strong ac-driving: a combination of Floquet 

and Van Vleck perturbation theory 
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We study the dissipative dynamics of a two-level system (TLS) exposed to strong ac driving. By 
combing Floquet theory with Van Vleck perturbation theory in the TLS tunneling matrix element, 
we diagonalize the time-dependent Hamiltonian and provide corrections to the renormalized Rabi 
frequency of the TLS, which are valid for both a biased and unbiased TLS and go beyond the known 
high-frequency and rotating-wave results. In order to mimic environmental influences on the TLS, 
we couple the system weakly to a thermal bath and solve analytically the corresponding Floquet- 
Bloch-Redfield master equation. We give a closed expression for the relaxation and dephasing rates 
of the TLS and discuss their behavior under variation of the driving amplitude. Further, we examine 
the robustness of coherent destruction of tunneling (CDT) and driving-induced tunneling oscillations 
(DITO). We show that also for a moderate driving frequency an almost complete suppression of 
tunneling can be achieved for short times and demonstrate the sensitiveness of DITO to a change 
of the external parameters. 

PACS numbers: 03.65.Yz, 42.50.Hz, 03.67.Lx 



I. INTRODUCTION 



The dissipative two-level system (TLS) has quite 
a rich and long history of both experimental and 
theoretical investigations [l[. Despite its simplicity, 
it is a very prominent candidate for modeling various 
different situations in physics as well as in chemistry 
and provides a testing ground for exploring dissipation 
and decoherence effects in genuine quantum-mechanical 
systems. The development of maser and laser technology 
triggered the examination of those systems under the 
influence of strong time-dependent driving fields, which 
yields dressed TLS states [2[ in turn leading to a variety 
of phenomena like coherent destruction of tunneling 
(CDT) |3H5| or driving-induced tunneling oscillations 
(DITO) For taking into account the influence of 

the environment, the driven spin-boson model 0, [T^] 
has proven to be a suitable candidate. 
In recent years, the driven TLS has experienced a strong 
revival in the field of quantum computation, as litho- 
graphic fabrication techniques allow the construction 
of artificial atoms that arc coupled to the modes of 
an oscillating field by a strength never reached in real 
atoms. Here the TLS implements the two logical states 
of a qubit. We mention just two prominent solid-state 
realizations of the qubit, namely, the Cooper-pair box 
[ill - [l4| and the Josephson flux qubit [l^l4l7j . Strong 
coupling between the TLS and a sing le oscillator photon 
was first successfully reported in [18| for a charge qubit. 
In this weak driving limit, the oscillator is usually de- 
scribed in its quantized version and occupied by a small 
number of photons Recently, a lot of theoretical 

effort has been put into solving the dynamics of such a 
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system [2Cj-|28|. Also for strong driving, the applied field 
can still be described by a quantized oscillator. However, 
for high photon numbers, the TLS-oscillator system is 
conveniently treated in the dressed state picture Q. In 
this strong driving regime, a series of experiments and 
theoretical investigations have been performed recently 
on superconducting qubits examining Rabi oscillations 
in the multiphoton regime and the validity of the dressed 
state picture d, [2fj| - |4G |. To account for environmental 
effects, the latter is usually combined with the phe- 
nomenological Bloch equations || H^, HH, |39| . 
With the applied field being in a coherent state and for 
high photon numbers, an equivalent description consists 
in replacing the quantized oscillator by an external, 
classical driving 0, |H|. Together with the coupling to 
a bath of harmonic oscillators, it leads to the driven 
spin-boson model [H, [i"oj . which has been examined by 
applying various techniques. For example, the (real- 
time) path-integral formalism provides a formal, exact 
generalized master equation for the dynamics of the 
reduced density matrix of the TLS, which can be solved 
approximately for certain parameter regimes nam- 
|48[. Among those treatments, the noninteracting blip 
approximation (NIBA) [l|, |4|| is the most prominent 
one and is based on an expansion to lowest order in 
the tunneling matrix element of the undisturbed TLS. 
It provides good approximate results for intermediate 
to high bath temperatures and/or strong damping of 
the system with arbitrary driving frequencies. However, 
at low temperature it fails to reproduce the dynamics 
of a biased TLS correctly. In (50rl53 |. the polaron 
transformation leads to an intcgro-diffcrcntial kinetic 
equation for the populations of the density matrix, which 
is equivalent to the generalized master equation under 
the NIBA. An alternative way to gain the dynamics 
of the driven spin-boson model for weak system-bath 
coupling and within the Markovian limit is to solve 
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the underlying Bloch-Rcdfield equations. This is done 
numerically for weak damping in 0, [H, [56| , while @, H(| 
additionally provide an analytical examination of the 
dynamics in the high-frequency regime. 
In this work we introduce a new approach to solving the 
dynamics of the monochromatically driven spin-boson 
model taking into account analytically the fast oscil- 
lations induced by the driving as well as the transient 
dynamics. In a first step, wc combine Floquet theory 
EE EH, H3 with Van Vleck perturbation theory @, [58| 
to derive the dynamics of the nondissipative system. 
This approach has recently been used also in [59| to 
evaluate the time-averaged transition probability of 
a nondissipative TLS. Going to second-order in the 
tunneling matrix element, we derive expressions which 
include the fast oscillatory behavior of the Floquet 
states and are beyond the common rotating-wavc results 
32 . 1 6 CM or perturbation theory in the driving strength 
4l| . |61| . Further, to analyze dissipative effects, we 



consider the regime of weak damping and solve the 
corresponding Floquet-Bloch-Rcdficld master equation 
applying a moderate rotating-wavc approximation. 
While in [62|, [63[ a similar approach is used to study the 
asymptotic dynamics of the driven spin-boson model 
perturbatively in the driving strength, our approach 
treats the full time evolution of the system, to all orders 
in the driving amplitude, in the regime of moderate as 
well as high external frequencies and for arbitrary static 
bias. Specifically, we are able to give closed analytic 
expressions for both the relaxation and dephasing rates. 
Our analysis enables us to shed light on the famous 
effects of CDT 0-@ and DITO Many investiga- 

tions of those phenomena have been performed in the 
high-driving regime. This work treats them analytically 
also for moderate driving frequency and amplitude. 
We examine both the nondissipative and dissipative 
cases and compare them to a numerical solution of the 
problem. 

The structure of the work is a follows. In Sec. |TT] 
the model Hamiltonian for the nondissipative system 
is introduced. We derive the corresponding Floquet 
Hamiltonian in Sec. Ill Al and analyze its quasienergy 
spectrum and the dynamics of the system in Sec. Ill Bl 
using a rotating-wavc approximation (RWA). In Sec. 
Ill CI we apply Van Vleck perturbation theory to second- 
order in the tunneling matrix element and compare the 
improved quasienergy spectrum to a numerical analysis. 
Further, we give in Sec. Ill Dl a detailed discussion of 
the parameter regime in which our approach is valid. 
To tackle the dissipative dynamics, we introduce the 



driven spin-boson Hamiltonian in Sec. Mil and solve 
the Floquct-Bloch-Redfield equation. We compare the 
analytical expressions for the relaxation and dephas- 
ing rates to the results obtained within the RWA and 
close the paragraph with a discussion of CDT and DITO. 

II. THE NONDISSIPATIVE SYSTEM 

In a first step, wc neglect environmental effects on the 
driven TLS and consider the Hamiltonian 



H T hs(t) = -- [Act x + (e + Acosujt)a z ] 



(1) 



Here, a z and a x are the Pauli matrices, and as basis states 
we choose the eigenstates of a z . | f) and | \) (localized 
basis) . The coupling strength A between those two basis 
states is time independent, whereas the bias point con- 
sists of the dc component e and a sinusoidal modulation 
of the amplitude A and frequency uj. 



A. Floquet Hamiltonian 

To resolve the dynamics of the driven system, wc take 
advanta ge of its periodicity and apply Floquet theory 
fiol [ill . l57j ] , about which we give a short overview in 
Appendix [XJ For the driven TLS it leads to the Floquet 
Hamiltonian "Htls (t) = #tls (t) — ihd t ■ Considering the 
case A = 0, we find the following set of eigenstates of 
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(2) 



with quasienergies hey^ n = =f|e — Hnu. Here, Jk{x) is 
the fcth-order Bessel function. In the composite Hilbert 
space T-L ®T [57( , which is introduced in Appendix [S] 
those states become 



K A ,n» = im>E J ±("-o (^)io, 



(3) 



with the state vectors { 1 1 ) } being a basis for T and | ai ) ) = 
| a) £§> \ l). For the case of a finite tunneling matrix element 
A, the Floquet Hamiltonian "Htls is nondiagonal in the 
above basis ([3]) and becomes in matrix representation 
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"Htls = ft 




(4) 



Wc defined 



the kind 



where we used the relation 164 



J n (u±v)= 2J Jn T k(u)Jk(v). 



A, 



(5) 



(6) 



To find the dynamics of the system, we have to diago- 
nalize the Floquet matrix. In the remaining subsections 
we discuss two approximation schemes. A rotating-wave 
approximation scheme is discussed in Sec. IIIB1 while in 
Sec. Ill CI Van Vleck perturbation theory is presented. 
We also show that the RWA results can be obtained with 
Van Vleck perturbation theory to lowest order in A. 



B. Rotating-wave approximation 

Let us look at the spectrum of the unperturbed prob- 
lem (A = 0). We notice that whenever the static bias 
fulfills the condition e = mui, the states \u° n )) and 
)) are degenerate, as then 



muj = 0. 



(7) 



In this case, we speak of an m-photon resonance. As 
long as A is only a small perturbation, u 3> A, then 
"Htls w iU exhibit a similar energy spectrum. The main 
corrections to the unperturbed Hamiltonian come from 
matrix elements connecting the (almost) degenerate lev- 
els. Thus, as a first approximation, we diagonalize an 
effective Hamiltonian, which consists of 2 x 2 blocks of 



-Ha 



-f A_ 



(8) 



and describes the energy states being connected by an 
m-photon resonance. 

This result is also obtained within the RWA scheme as in- 
troduced in {32l 60]. In those works, the time- dependent 
system Hamiltonian (JTJ is transformed to a rotating 
frame, and only terms fulfilling the resonance condition 
(|7|) arc kept, while the fast-rotating components are ne- 
glected. This RWA is different from the conventional 
Rabi rotating-wave approximation, which is perturbative 
in the driving amplitude A, see, e.g., @,Hl|, and becomes 
exact for circularly polarized radiation. In contrast, the 
RWA we are using treats the driving amplitude nonper- 
turbatively. 

Concerning the eigenenergies of the Floquet Hamilto- 
nian for finite A, we notice that the exact crossing of 
the unperturbed energies (A = 0) at e = mui becomes 
an avoided crossing (see Fig. [1} and the perturbed eigen- 
states are a mixture of the unperturbed ones. Those with 
higher eigenenergies are labeled |3>+^+ m )); those with 
lower energies |$?:^ A )) They are defined below. 

In the far off-resonant case, |$5^1 m )) corresponds for 
e > muj to the unperturbed state |ti? „ +m )), and \$™J^)) 
to \uS. n )). For e < mui, the state |$^n+m)) corresponds 
to and |$^ A )) to |u° i „ +m }}. The eigenenergies 

are 



RWA 



he RWA 



(9) 
(10) 
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For the special case of vanishing static bias (e = ) and 
0-photon resonance, 



ifJ?f(*) = cos- 8 [\J (A/u;)A\- 



t 



(16) 



which agrees with the high-frequency result, u 3> A, of 
earlier works SO El. 



C. Van Vleck perturbation theory 



FIG. 1: (Color online) Quasienergies and e+^ A m (trian- 

gles) and unperturbed quasienergies and s^ n+m (dashed 
line and dotted-dashed line) for an m-photon resonance. The 
unperturbed quasienergies show an exact crossing at e = mu) 
according to Eq. @. The corresponding eigenstates are 
|-it° n )) and \u° tn+m )). For finite A an avoided crossing can 
be observed. The energy e5^ A and the corresponding eigen- 
state |<1? I !^ A )) are represented by black downward triangles, 
whereas e+^ A m and |<E>+^ A m )) are shown by red upward 
triangles. For e > mui we find that |$?:^ A )) approaches 
|w?,n))) while |$lJ^ A _ m )) becomes |uJ |T i+ m )) and vice versa 
for e < mu). The labeling of the perturbed eigenstates and 
eigenenergies is chosen in a way that £+^+ m > e-^ii for all 
e. 



with the oscillation frequency 



(11) 



The corresponding eigenstates are 
|^ A »=-sin^|4„)) 



I^J)=co S ^-K, n » 



qRWA 

sign (A_ m ) cos ^_|4„ +m )), (12) 

qRWA 



RWA 
*n L.O 



sign (A_ m ) sin -^_|^ n+ J), (13) 



where 



tan9* WA 



|A_ 



for 0<e^ A <7r. (14) 



-e + muj 



Now we are able to recover the time-dependent dynamics 
of the system (see Appendix [C]). As an example, we give 
the state returning probability for a system starting in 
the localized state | I) and returning to this state: 



P^ A (i) = cos 2 ( 



cos 2 65; WA sin 2 a RWA 



t 
2 
(15) 



As pointed out already in [32], (59[, the RWA fails in 
explaining higher order effects in A such as a shift in the 
oscillation frequency. Furthermore, we will show that 
the couplings between the nondegencrate states in (TJ| 
are needed to get physically correct expressions for the 
relaxation and dephasing rates. In the following, we will 
use Van Vleck perturbation theory to go beyond those 
shortcomings. Originally this method was used to treat 
modifications on diatomic molecules caused by vibrations 
and rotations of the nuclei (H)]. Since then the formal- 
ism has found many applications in both chemistry and 
physics and experienced various modifications; see, for 
example, [65M68| . The main formalism behind these dif- 
ferent variants is, however, always the same: a unitary 
transformation U is applied in order to construct an ef- 
fective Hamiltonian which exhibits, to a certain order in 
the perturbation, the same eigenenergies as the original 
Hamiltonian but only connects almost degenerate levels. 
In this work we choose for the transformation the form 
U = cxp(iS'), which was originally proposed by Kem- 
ble in [69| and is described in more detail in 0. In the 
case of the Floquct Hamiltonian, the effective Hamilto- 
nian then becomes: H c s = cxp(iS')'HTLS exp(— iS). Wc 
calculate the transformation matrix S up to second-order 
in A. As shown in Appendix iBl the so-obtained effective 
Hamiltonian for an m-photon resonance again consists 
of 2 x 2 blocks, however, compared to the one of the 
previous section, it has corrected diagonal entries: 
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It leads to the new quasienergies 



(17) 



te_ in = ft[(-n-im)w-inW] ) (i 8 ) 
fie+, n+m = H[(-n - i m )w + ^nW], (19) 



with the second-order oscillation frequency [81 



\ 



2 , e + Iuj 

l^t — m 



A 2 



(20) 



Compared to the frequency obtained within the RWA, 
Eq. (jTTJ) , this new frequency is shifted due to the second- 
order elements in (|17p , and the condition for an m-photon 
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resonance reads now 



S = Did) 



1,5 



2 , j-^ e + luj 

ljt—rn 



(21) 



In Fig. [3J we compare Eqs. (JTSJl and (fTi?j) for the 
3 




FIG. 2: (Color online) Quasienergies against static bias e. 
The upward triangles result from numerical diagonalization 
of the Floquet matrix Q , while the solid lines correspond to 
the analytical formulas (|18[) and (|19[) . Parameters are cj/A = 
2, A/A = 3. At an m-photon resonance, we find avoided 
crossings with a gap distance of A m . 

quasienergies against the eigenenergies we find from nu- 
merical diagonalization of the Floquet matrix (Q| . When- 
ever the resonance condition [Eq. ([21]) ] is fulfilled, we 
notice avoided crossings whose gap distance is deter- 
mined by A m for an m-photon resonance. The eigen- 
states |$± n)) of the effective Hamiltonian are the same 
as in ([T2"j) and (fT3")l , with the mixing angle 6™ replaced 
by 



e,. 



arctan 



2 



— m s+Ilu 



To get the eigenstates of "Htls, we calculate |$±. n )) = 
cxp(— | <I>5? Tl ) ) and, following Appcndix[C] the survival 
probability P±_y±(t). In Fig. [3l we visualize the results 
for the survival probability close to a 2-photon resonance 
obtained from the RWA approach and first- and second- 
order Van Vleck perturbation theory. We notice that by 
applying the RWA the fast oscillations in the first-order 
Van Vleck result are averaged out. When we compare 
first- and second-order predictions, the shift of the oscil- 
lation frequency is striking. But also the amplitude of the 
oscillations changes, which is due to the corrected mix- 
ing angle, Eq. (|22j) . Inserting the second-order mixing 
angle and frequency into the RWA formula (| 15[) results 
in averaging over the fast oscillations of the second-order 
Van Vleck graph. To also cover the fast driving-induced 
oscillations, it is essential to use the eigenstates |3>± n )) 
instead of the effective ones, which leads to a more com- 
plicated formula for P^j_(t), sec Appendix [Cl 



■-■ RWA 

1st order VV 

averaged 2nd order VV 

2nd order VV 




FIG. 3: (Color online) Survival probability Pj__>|(t) close 
to a 2-photon resonance. The parameters are e/A = 4.1, 
cj/A = 2, A/A = 3. We compare results obtained from a 
RWA and the first- and second-order Van Vleck perturbation 
theory. Furthermore, we show the averaged second-order Van 
Vleck dynamics, which correspond to Eq. (|15|l with the Van 
Vleck frequency [Eq. (|20|l ] and second-order mixing angle 
[Eq. I®]. 



D. Validity of the Van Vleck approach 

In closing this section, we give an overview of the pa- 
rameter regime in which our approach is valid. To ap- 
ply Van Vleck perturbation theory at all, a requirement 
for the Floquet Hamiltonian is that it has for finite A 
a similar doublet structure as in the unperturbed case 
(A = 0). This means that the off-diagonal elements in 
([4]) connecting different doublets with each other must be 
much smaller than the distance between those doublets 



2]: 



m+n/ 



^ \ e t,n ~ £ l,l+m+n\ 



(23) 



(22) for any 1^0. Using Eqs. ([5]) and ((TJ), this becomes 



|A_,_ r 



< £ 



(J 



(24) 



Because |A_;_ m | < A, this condition can be even ful- 
filled for frequencies lu < A. Once Eq. ([24]) is valid, we 
still have to check at which order one can stop the pertur- 
bative expansion in A. We will distinguish now between 
two situations: the case of being close to or at an m- 
photon resonance and the regime far from resonance. 



1. Dynamics close to or at resonance 
Using e ~ mw, Eq. ([24]) becomes simply 



|A_ 



■l—m | 



(25) 



Notice that the right-hand side of (f2~5]) still depends on 
Ajuj. Thus, while being surely fulfilled in the RWA case, 
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ui ^> A, condition ([25]) is in general less restrictive. To 
show this, we examine the following two limiting cases. 
First, the limit A/us 1 is considered. For arguments 
with < x <C \fn + 1 .the nth-order Bessel function 
becomes approximately [7 



\Mx)\ 



2l"l|? 



Thus, for Aju) <C 1, we find that 



and (|25p becomes 



1 2H|n|! 

(A/w)l-'- m l 



w > A 



(26) 



(27) 



(28) 



Because A/w < 1, Eq. ([28]) is fulfilled for any Z ^ if it 
is satisfied for I = —in: i.e., if 



A 

Iml 



(29) 



In the case of a 1-photon resonance, this leaves us with 
the RWA condition, uj A, as then nearest neighbor 
doublets are connected by a Ao element in the Floquct 
matrix which approaches A for small A/uj. All other per- 
turbative off-diagonal entries in (U]) are vanishingly small. 
In the case of an m-photon resonance with m ^ ±1, the 
dressed element A connects more distant doublets, so 
that the Van Vleck condition ([25]) can be realized ac- 
cording to (|2"9"|) for frequencies smaller than the ones de- 
manded by the RWA. 

In the opposite limit of A/ to ^> 1, an upper bound for 
the dressed Bessel function is [7fl 



|A„| < Av^VA. 
Using this, we find that (j2"5j) is verified if 



A2 
A 



(30) 



(31) 



Since A > w, it follows that A 2 <C Alo < A 2 and thus 
A ^> A. Hence, Eq. (f3"Tj) represents an improvement to 
the RWA condition. 

Further, being close to an m-photon resonance, one sin- 
gle frequency will dominate the system's behavior, and 
thus neglecting the remaining fast-oscillating terms will 
already give a good picture of the coarse-grained dynam- 
ics. This dominating frequency is represented by fi^ WA 

(2) 

and by flm in the case of the RWA and the second-order 
Van Vleck perturbation theory, respectively. To obtain 
those frequencies, it is enough to diagonalize the corre- 
sponding effective Hamiltonian, without yet considering 
any modification of the eigenstates of the effective Hamil- 
tonian. As shown in the previous subsection, ft* WA cor- 
responds to the main frequency of the system obtained 



by applying Van Vleck perturbation theory to first order 
in A. Naturally the question arises as to how good these 
approximations are, or which orders in A are necessary 
depending on the parameter regime. 
In a first step, we examine the improvement obtained by 
using second-order Van Vleck perturbation theory com- 
pared to the RWA; that is, we consider 



0.15 




0.05 



FIG. 4: (Color online) Comparison of the main oscillation 
frequency obtained by the RWA and second-order Van 

(2) 

Vleck frequency Hi, for a fixed static bias, e = 4. OA. The rel- 
ative mistake performing the RWA is shown against the driv- 
ing frequency cj/A and driving amplitude A/ A. The darkest 
areas show regions in parameter space of small or no devia- 
tions between the two approaches, whereas the lightest areas 
show a deviation of 15 % or more. 



,EWA 



(2) 



(2) 



(32) 



and plot it in Fig. 0] against the driving frequency u> and 
amplitude A at a fixed value of the static bias e. The de- 
viations are visualized through different shades of color. 
The lightest areas stand for a relative mistake of 15 % or 
more. We can tell from Fig. @]that the RWA fails for low 
driving frequencies and/or weak driving amplitudes. The 
darkest areas determine regions in the parameter space 
where almost no difference between the RWA and second- 
order Van Vleck approach can be found. Of course this is 
no indication that the results are reliable in those areas, 
but rather that second-order perturbation theory yields 
no improvement to the RWA. 
To check the accuracy of the second-order Van Vleck 

(2) 

frequency Qm , we calculate the deviation 



,(2) 



l.o 



(2) 



(3) 
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(33) 
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g -| 5 (t) an d its Fourier transform 



0.05 



10 15 20 25 

A/A 

FIG. 5: (Color online) Comparison of the main oscillation 
frequency Qm obtained by second-order and third-order Van 
Vleck perturbation theory for a fixed static bias, e = 4. OA. 
The relative mistake performing second-order perturbation 
theory is shown against the driving frequency u/A and driv- 
ing amplitude A/ A. Color scale is the same as in Fig. [4] 



F{y) 



(34) 



at resonance e — uj but for a driving amplitude with 
|Ji(A/w)| 7^ 0. One clearly sees that one frequency, 
namely, f2i, is dominating, and already the RWA con- 
veys a good impression of the dynamics. 




FIG. 6: (Color online) Survival probability P±_>±(t) for e/A = 
4, oj/A — 4, and A/ A = 4.1. Exact numerical results are 
compared with RWA and second-order Van Vleck results. 



(3) 

from the frequency VL m obtained applying Van Vleck 
perturbation theory to third order [5£ 



Results for e^ 2 ' 
are shown in Fig. [5] Again we only consider mistakes 
up to 15 %. We find strong deviations in the region 
of low driving frequency and intermediate driving am- 
plitudes. In the remaining parameter space, the agree- 
ment between second- and third-order Van Vleck per- 
turbation theory is quite good apart from small islands. 
Those islands are located at values of uj and A where the 
second-order condition for coherent destruction of tun- 
neling (CDT) is fulfilled, see discussion in Sec. IIII CI 
For example, for e/A = 4.0 and ui/A = 2.0, they oc- 
cur at the zeros of the Bessel function J^A/ut). Since 

(2) 

at those points the second-order frequency £l m vanishes, 
even small third-order contributions yield a significant 
correction. This behavior visualizes nicely the findings 
of Barata et al. [7l| and Frasca [7^, who proved ana- 
lytically that fl m does not completely vanish at the ze- 
ros of the Bcsscl function if third-order contributions in 
A are taken into account. On the contrary, both the 
RWA and second-order Van Vleck perturbation theory 
predict a vanishing frequency at those points and there- 
fore agree perfectly with each other in Fig. |U We want 
to emphasize again that, as can be seen from Fig. [5j our 
approach also yields good results for low driving frequen- 
cies, uj < A, and small driving amplitudes, A ~ A. 
In Figs. [5] and [7J we show the survival probability 
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FIG. 7: (Color online) Absolute value of the Fourier transform 
F(u) of Pi^i(t) in Fig. [6] The oscillation frequency corre- 
sponding to fii = 0.45A is dominating and is also predicted 
by the RWA approach. Notice, that there is exact agreement 
between the second-order Van Vleck and the numerical re- 
sults. 



2. Dynamics away from resonance 

The situation changes when we are away from a res- 
onance. Already intuitively it becomes clear that the 
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dynamics will not be governed anymore by a single fre- 
quency. Therefore, by looking only at the coarse-grained 
dynamics of the system and averaging out the driving- 
induced oscillations, significant information is lost. This 
case is presented in Figs. [8] and HI where we are in the 
region between the 1- and 2-photon resonances. In con- 
trast to Figs. |5]and [3 we find that several frequencies are 
dominating and determine the dynamics of the system. 
The second-order Van Vleck approach reflects this be- 
havior almost perfectly because the driving-induced os- 
cillations are accounted for. However, the RWA shows 
only one single oscillation because the others are aver- 
aged out. It depends on the choice of m in the formula 
for the RWA, Eq. (fl~T|) . which of the frequencies is taken. 
This explains also the cuts in Fig. |4] see, for example, 
the horizontal line just below uj/A w 3. At these values 
of the frequency, we change m in our analytical calcu- 
lation. In Fig. [5j those cuts are barely visible. Being 
away from the resonance point, the modifications of the 
external driving on the system's eigenstates must not be 
neglected. 

Off resonance, the requirement (|24[) for Van Vleck per- 
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FIG. 9: (Color online) Absolute value of the Fourier transform 
of Pi->±(t) in Fig. [8] The numerical and the second-order 
Van Vleck graphs clearly show several dominating oscillation 
frequencies, while the RWA only shows one of them. Note 
also the different scale of the y axis compared to Fig. [7] 




FIG. 8: (Color online) Survival probability P|_->j, for e/A = 4, 
ui/A = 2.7, and A/A — 4.1. Exact numerical results are 
compared with RWA and second-order Van Vleck predictions. 



turbation theory is surely fulfilled for a large enough 
static bias, 



IA. 



< £. 



(35) 



X]fc fcvkoLbki which are bilinearly coupled to the TLS by 
the coupling Hamiltonian 



Hi, 



k 



(37) 



III. THE DISSIPATIVE SYSTEM 



Here x = u z /2 is the position matrix of the TLS and 
Afc the coupling strength to the fcth mode of the bath. 
The spectral density of the bath can be expressed as 
G{ v ) = Ylk^k( y ~ v k)- We assume further that at 
time t = the bath is in thermal equilibrium and un- 
corrclated to the system, so that the full density matrix 
W(t) associated with H(t) has at initial time the form 
W(0) = p(0) (g> pb(0), where pit) is the density matrix of 
the TLS and p B (0) = cxp(-/3ff B )/tr B cxp(-/3 J ff B ) is the 
density matrix of the bath at temperature T = [ksP)^ 1 . 
Following [z34z3 and performing a Born and Markov 
approximation, we arrive at the Floquet-Bloch-Redfield 
master equation 

Pap(t) = -i( £ a - ep)p a p{t) + 7T ^2 £aP,a'0'{t)p a >p>(t), 

a',P' 

(38) 

where the density matrix is expressed in the basis of the 
energy eigenstates of the TLS: 



To include dissipative effects on our system, we con- 
sider the time-dependent spin-boson Hamiltonian [l], H, 


H(t)=H TLS (t) + H B + H int , (36) 

where the environmental degrees of freedom are mo- 
deled by an infinite set of harmonic oscillators, H B = 



Pa p{t) = {§ a (t)\p(t)\$f,{t)), a,0 = ±. (39) 

Notice that e a = e a fi and \<fr a {t)) = |4> a>u (t)). Correc- 
tions to the oscillation frequencies due to the Lamb shift 
arc not accounted for. The first part of (|3"ci)) describes 
the nondissipativc dynamics as treated in Sec. [TTJ The 
influence of the bath is fully characterized by the time- 
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dependent rate coefficients 



e i("+" MJ (Naa',n + Npp' ,- n ')X aa > , n X p' p y 

— Spp' X a p" <n i Npiigi , n Xp'> a i t n 
P" 

^aa' ^ "j a" j3' n-^ j3' a" ,n-^- a." f3 ,7i' r 



(40) 



with N a p <n = N(e a - sp + nuj), N(is) = G(v)n th (v), 
and n t h(i>) = \ [coth(/i/3^/2) — 1]. As also the ma- 
trix elements of the position operator x = a z /2 are 
periodic in time, we express them in a Fourier series, 
{$ a (t)\x\*f,(t)) =E„e i "-*X Q/3 ,„. 



A. Position matrix elements 

The Fourier coefficients appearing in the rate equations 
40| can be calculated by 

f 



X a p tl = -J dt e -""(* a (t)\x\*p{t)) 
= (($ a ,oM$ i 8 I i», 



(41) 



where we used the periodicity of the eigenfunctions of 
the TLS and the definition of the internal product in the 
composite Hilbcrt space, Eq. (|A5|) . From this we find 
that X a p^ n = Xp a n and that we can use the Floquct 
eigenstates (|B12j) and (|B13|) to calculate the Fourier co- 
efficients to second-order in A. We get 



— hn(®m 

sine 



E 



k—m 



A_ 



£ + (n — k — m)u>] [— e + kto] ' 

(42) 



X 



(2) 



X^JQ, 



cos 6 I; 
8~~ 



E 



A„_A.A_ t 



e + (n — k)uj] [— e + kuj] 



(43) 



with 



(i) _ sin£ sign(A_ m ) 



. 2 £ A_„ 
sin - - 



2 £ A n _2m 
COS - - 



2 e + (n — 2m)uj 



2 — e + nui 
(l-8n, m ), (44) 



Y (i) f n cos£ sign(A^ m ) . 

= ^— On,o + 1 sin£ 

A_T7i_77 An— m 



-e + (to + n)cj e + (n — to)w 



(1 - <*n,o) 5 
(45) 



where for £ cither the mixing angle 6^j WA or m is 



used. Further, we find that X 



(2) 

++,» 



-Xi 2 l„. Within 



the RWA, we would get A^a = i sm eRWA^ and 
= ~\ cose^ WA (5„ i0 - From this we notice that, in 
the case of a simple RWA, A^+A would be nonzero for 
n = to and X^f^ for n = only. An improvement to 
that can already be achieved by using Van Vleck pertur- 
bation theory to first order in A, yielding x£p „(6™). 
It contains next to the RWA results additionally first- 
order corrections for any index n in X a p iU . Figures [TOl 
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FIG. 10: (Color online) Fourier coefficient \X >n | for various 

values of n against driving amplitude A. We examine the case 
e = 2cj. The parameters are e/A = 4.0 and uj/A = 2.0. The 
black triangles show data points from numerical diagonaliza- 
tion of the Floquet matrix; the red (dark gray) dashed curve 
is obtained from first-order perturbation in A [Eq. ()45f) ] . 
whereas the green (light gray) solid curve is obtained by going 
to second-order in A [Eq. (|43[) ]. 

and ITTI show the absolute value of the coefficients X 

and X h) „, respectively. We find a good agreement be- 
tween the results obtained by a numerical diagonalization 
of the Floquet matrix ([4]) and second-order Van Vleck 
perturbation theory, Eqs. ([4"2"| and (|4"3"| . Concerning 
Figs. UOf b) and (d) we see a qualitative improvement 
by going from first to second-order in A. While in Fig. 
UOf b) the first-order result approaches a nonvanishing 

coefficient X 2 f° r A — > 0, Eq. (|4"3"| corresponds to 

the numerical calculation very well even in the region of 
low driving amplitude and meets our expectation that 
all Fourier coefficients except for X a p,Q vanish at zero 
driving. The problem of the first-order results at low 
driving strength is caused by the definition of the first- 
order mixing angle 9^j WA , Eq. (fl"4|) , which is 7r/2 for 
e = muj. When n ^ ±m in A a l „(6^ WA ), the coeffi- 
cient approaches zero for A — > because of the term 



A_ r 



— £ + (to + n)uj £ + (n — m)uj 



(46) 



in Eq. (|45j) . However, for n = ±m a zeroth-order Bessel 
function occurs in that part which does not vanish for 
A — > 0. A second-order improvement of the mixing angle 
as done in Eq. (j2"2"j) solves this problem. In Fig. HUf d) the 
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FIG. 11: (Color online) Fourier coefficient \X |_ jTl | for various 

values of n against driving amplitude A. The parameters are 
the same as in Fig. 1101 



first-order solution predicts a coefficient xi 1 ! o(0^i WA ) 
which is constantly zero. 

Also, in Fig. [IT] a noticeable improvement between first- 
and second-order perturbation theory can be seen. In 
ITTT c) the first-order solution shows a constant coefficient 

X_1 2 (G^ VA ) = 0.5. We see from the numerics and 
second-order results that indeed the constant value is 
reached asymptotically for high driving amplitudes; how- 
ever, for small driving amplitudes, we find a vanishing 
coefficient. In Fig. [TTTh) we can observe a behavior like 

infTOTb). namely, that X^\_ 4 (8jj, WA ) does not approach 
zero for A — > 0. The explanation is similar to the above 
case. 



C 



MKWA 



v. 



a.' ,n-X 



8" a.' ,n 



a." ,nX a " } _ n / . 



(47) 



We observe that £^ Q R ^ 



/■MRWA ^MRWA _ fMRWA 
^aa.-\ — ? — \-,aa — ,flai 

/-MRWA _ /'MKWA nrlr | /-MRWA _ /-MRWA Mnrprwpr 

P— (f) + p++(t) = 1 and p+_ (t) = pl+(i). This yields 
simple expressions for the reduced density matrix ele- 
ments to first order in the coupling n to the bath: 

£MRWA 

p__(t) = 7T ~~' ++ + ^-(mw + fiMje-'W 

7rcl 7T 

+ 2£ MR ^A Rejcdephe-^+^^le-^ 4 , (48) 



p_ + (i) = Crol (£™ - £MRWA_ )e - 7rolt 

1 ^MRWA , jMtp ->H+tlL 2) )* p -7d, 



T 2^ — |-.H — ^dephe 

+ c d cph i (^ + ^y^+^^e-^ 4 . (49) 

7T 

The constants c re i and Cdoph are fully determined by the 
initial conditions, Eq. (|C9p . The expressions for the 
relaxation and dephasing rates are 

//.MRWA cMRWA \ /.MRWA 
7rel = 7T (1, h+ — i, j , 7deph = — 7Ti-__| ,. 

' . (50) 

With (|47[) we can express them in terms of the position 
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B. Moderate rotating- wave approximation 

Having calculated the position matrix elements, our 
rate coefficients C a /3, a '/3'{t) are fully determined. What 
remains to do is to solve the Floquet-Bloch-Redfield mas- 
ter equation (|38p for the density matrix p. For an an- 
alytical calculation, there is, however, still a difficulty: 
the time dependence of the coefficients. To get rid of 
this, we perform a moderate rotating-wave approxima- 
tion (MRWA) [75| : i.e., we neglect fast-oscillatory terms 
in (I40p . which amounts to selecting only the terms with 



FIG. 12: (Color online) Relaxation (a) and dephasing rate 
(b) against driving amplitude A for u/A = 2.0, e/A = 4.1, 
h/3A = 10, and k = 0.01. Results obtained within the second- 
order Van Vleck perturbation theory are compared with RWA 
calculations. Notice that the RWA predicts an unphysical 
vanishing of the relaxation rates at the zeros of A_2- 



matrix elements, yielding 



7rel 



4-E 



X: 



(51) 
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7doph = ^7rcl + 4?r^ N — , n X 2 __ > 



(52) 



With (|42|) and (|43|) we arrive finally at one major result: 

7rcl = 7rcl + 7 "ol and 7deph = 7deph + 7 dcph 

(53) 



with the contributions 



7rel = 



\ coth (MfiO 



'-si! 



A 2 



(54) 



7 ; l cl = ttG 



coth 



h/3 
~~2 



— sin 



2 ®m ^ — (n+m) 

2 e—(?i + m)w 



COS 



2 ®m 



2 e + (n — m)cj 

(55) 
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FIG. 13: (Color online) Survival probability Pj,_ > j_(i) close to 
a 2-photon resonance. The parameters are e/A = 4.1, cj/A = 
2.0, A/A = 3.0, ft = 0.01, and ft/3A = 10. Analytical results 
obtained by second-order Van Vleck perturbation theory are 
compared with RWA results. The inset shows the long-time 
dynamics and visualizes the deviation of the RWA from the 
Van Vleck dynamics in the asymptotic limit. 



and 



7dcph = ^7rcl + TTN{0) COS 2 Q r 



1 - 



A 2 



2 (e - few) 2 



(56) 



"Tdcph 



r G(nw) 



hp 

com — nui 
2 



- 1 
A„ 



sin 2 e, 



-e + (m + n.)w — e + (n — m)ui 



(57) 



For zero temperature, an instructive interpretation of 
those rates in terms of a dressed energy level diagram 
is given in (39j. Within the RWA, on the contrary, the 
corresponding rates read 



RWA 
7rel 



nG coth sin 2 0™* (58) 



and 



RWA 1 „,RWA , ^ w /n\ „ 2 oRWA /r \ 
7deph = ^rel + ttA (0) cos O m . (59) 

The RWA rates correspond to those of an undriven TLS 
[H using the dressed energy levels and the RWA mixing 
angle 6* WA 

In Fig. 1121 we compare the rates obtained through Van 
Vleck perturbation theory with the RWA ones for an 
Ohmic spectral density, Gohm(^) = Kl/ , where K is the 
dimensionless coupling constant between TLS and bath. 
For both the relaxation rate, Fig. [T27ab and the dephas- 
ing rate, Fig. [T2T bb the RWA approach underestimates 
the rates. The failing of the RWA becomes especially 
evident in Fig. I12f ab where a zero relaxation rate is 



predicted for driving amplitudes A under which A_ m 
vanishes. This implies in particular no relaxation at zero 
driving and m 7^ 0. Again we see that the higher order 
matrix elements in the Floquct matrix Q are necessary 
in order to correctly describe the dynamics. We find 
that for certain values of the driving amplitude, namely, 
whenever A_ m = 0, 7^ vanishes and thus 7 ro i becomes 
minimal, a behavior which could be already predicted by 
inspecting formula (|42|) for A_ +i „. This could be ex- 
ploited experimentally to minimize relaxation. On the 
other hand, for higher driving amplitudes, 7d G ph exhibits 
peaks at A_ m = because of the cosine in j^ oph - For 
a high driving amplitude, our rates approach asymptot- 
ically the ones predicted by an RWA approach. In the 
opposite regime of small driving amplitudes, however, de- 
viations between the RWA and Van Vleck rates occur, as 
matrix elements connecting the different doublets in the 
Floquct matrix play a more important role. Common to 
both approaches is that the external driving yields a re- 
duction of the rates with increasing strength, a behavior 
which was already numerically predicted, e.g., in fbH ]. 
The failure of the RWA also becomes evident in Fig. [TBI 
where we show the dissipative dynamics obtained for an 
Ohmic environment. Comparing the results for Pj,_ > |(t) 
which wc obtain from second-order Van Vleck perturba- 
tion theory - formulas (|4"8"1) and (|4"9"|) combined with (JC2J) 
- with the RWA result, we find striking differences. Con- 
sidering the long-time dynamics (the inset in Fig. I13|) . 
we see that the RWA predicts quite a different asymp- 
totic value for Pj-^t). We notice further that the RWA 
exhibits a single oscillation frequency, which decays com- 
pletely to a constant value, while within the Van Vleck 
solution P|_^(t) oscillates for t — s- 00 around the equi- 
librium value. This latter behavior corresponds to the 
continuous driving of the system through the external 
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FIG. 14: (Color online) Absolute value of the Fourier trans- 
form F(u) of the survival probability for the RWA and second- 
order Van Vleck perturbation theory. The parameters are 
the same as in Fig. 1131 Next to the relaxation peak at 
v = the RWA dynamics are governed by a single frequency 
v = Q 1 ^ ffA . The second-order dynamics also exhibit the relax- 
ation peak and a main frequency, which, however, is shifted 

(2) 

to v = fim . Additionally, the higher harmonics of the driving 
can be seen in the second-order dynamics. For visualization of 
the S peaks, appearing at v = nu), a finite width and height 
have been artificially introduced. Furthermore, broadened 
peaks appear at v = nuj ± fi^ ■ 



field. It is completely missed by the RWA approach. 
For a further analysis of the dynamics, it is helpful to 
consider the Fourier transform of P^^it), see Fig. [14] 
Both the RWA and Van Vleck dynamics exhibit a relax- 
ation peak at v = and the dressed frequency of the 
system at v = S1„ WA and v = Qm , respectively. Those 
latter peaks have a finite width due to the dcphasing. 
Within the RWA f2^ WA is the only frequency; while for 
second-order Van Vleck dynamics, we find additional fre- 
quencies. They result from the higher harmonics of the 

driving and are located at integer multiples of the driv- 
el 

ing frequency, v — nuj, and at v = ncu ± Q m . The peaks 
at v = nuj arc S shaped as they suffer no dephasing, 

(2) 

whereas the peaks at v = nu ± fim show the broaden- 
ing of the main frequency. Already in the nondissipative 
dynamics, Eq. (|C17|) . we found the appearance of those 
multiple frequencies. They result from the beyond-RWA 
contributions in (|C12[) and (|C13|) and reflect the external 
driving. Dephasing only influences the dressed frequency 

£1$ in p (t) and P-+(t), see (gg} and (gSJ), and thus 

for the equilibrium state, the laser frequency at v = nuj is 
dominating. This asymptotic behavior agrees well with 
the findings in [H, S3, 



C. Coherent destruction of tunneling 

It has been found in 0, H| for a driven double- well po- 
tential and for a driven TLS in Q that under certain con- 



ditions, coherent destruction of tunneling (CDT) occurs. 
For a symmetric TLS (e = 0) and for high enough driv- 
ing frequencies (a; 3> A) this phenomenon was predicted 
to happen approximately at the zeros of Jq(A/lj), as can 
also be seen from Eq. (fT6|) . For a nonzero static bias 



and high frequencies, the necessary conditions for CDT 
are e = mw and J m (A/uj) = [7g, |77|. In this section, 
we compare the predictions of the RWA and Van Vleck 
perturbation theory against exact numerical results. 
For the case of an exact m-photon resonance (e = 
moj) and nonvanishing A_ m , the RWA mixing angle is 
e^ WA = tt/2, and we get from Eq. (fI5]t. 



J- m {A/u)A\ 



(60) 



Also from this formula, we see that CDT occurs at the ze- 
ros of J- m (A/uj). Notice, however, that for J_ m (A/uj) = 
0, Eq. dT5j) predicts P™^(t) = 1 even for systems which 
are not at an m-photon resonance; i.e., within the RWA, 
the condition e = muj is not necessary for CDT. 
Interestingly, also second-order Van Vleck perturbation 

(2) 

theory predicts Qm = for e = muj and J- m (A/uj) = 0, 
see Eq. ([20]) . However, as shown in [7l|,[72| and discussed 
in Sec. Ill Dl this condition holds only to second-order in 

(3) 

A; third-order corrections will cause ttm to be small but 
finite for e = muj and J- m (A/u) = 0. Thus, instead of 
being localized, the dynamics will oscillate with a large 
period. 

To visualize this behavior, we examine in the following 
without loss of generality the case of a 3-photon reso- 
nance. We choose w/A = 2.0 and e/A = 6.0. Then the 
first zero of J_ 3 (A/tu) occurs at Aj A w 12.7603. Using 
those parameters in Eq. (|20l) . the Van Vleck oscillation 

frequency Q 3 is zero. Figure [15] shows a comparison 
between the RWA and Van Vleck dynamics to second- 
order and an exact numerical treatment of the Floquct 
Hamiltonian for the above parameters. For the RWA, 
we see a complete destruction of tunneling because the 
driving-induced oscillations are not accounted for. Also, 
within the Van Vleck description, the coherent oscilla- 
tions are strongly suppressed; however, we notice fast 
oscillations because of the external driving. This be- 
comes especially clear in Fig. IToTb). At t = (2n + 1)tt/lu 
with n = 0,1,2,3..., we find sharp dips. The plateaus in 
between show weak oscillations whose number changes 
with m. The situation changes strongly for the numeri- 
cal graph: instead of a localization, a complete inversion 
of the population occurs; CDT seems to have vanished 
completely, as VL m is not vanishing. Considering, how- 
ever, the time scale in Fig. ITST a). we notice that the pe- 
riod of 27r/f2 m is rather large. For short times, see figure 
I15f b). also the numerical dynamics appear to be local- 
ized. Note that this observation also holds for the high- 
frequency case examined in Q : considering the dynamics 
at long times the localization will also be destroyed there 
due to higher order effects. 

In Fig. 1161 CDT under the influence of dissipation is 
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FIG. 15: (Color online) Coherent destruction of tunneling for the nondissipative case (ft = 0). The survival probability Pj__s.^(t) 
is shown at a 3-photon resonance. The parameters are e/A = 6.0, ui/ A = 2.0, and A/A = 12.7603. The Van Vleck solution is 
compared with the RWA and a numerical diagonalization of the Floquet Hamiltonian. Within the RWA, a complete destruction 
of tunneling can be observed, whereas the analytic Van Vleck solution exhibits driving-induced oscillations. The numerical 
solution predicts, on the contrary, complete population inversion with low but nonvanishing frequency Q m . Figure (b) is a 
blowup from figure (a) for a shorter time scale. There, the numerical and Van Vleck solutions agree well, and one can nicely 
see the small oscillations resulting from a 3-photon absorption or emission. 
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FIG. 16: (Color online) Coherent destruction of tunneling for the dissipative case (k = 0.01, h/3A = 10). The remaining 
parameters are the same as in Fig. [T5J A comparison between the numerical solution of the Floquet Hamiltonian and the full 
master equation (|38|l . the Van Vleck combined with the MRWA approach, and RWA results is provided. The RWA approach 
predicts a slower relaxation than the numerical one and Van Vleck solution. Deviations between the numerical and Van 
Vleck solutions can be seen especially in the long time limit. As in Fig. 1151 the numerical result predicts oscillations with a 
nonvanishing frequency f2 m . For short times, see figure (b), the numerical and Van Vleck results agree well. 



examined. As in Fig. 1151 we investigate a 3-photon res- 
onance with vanishing frequencies f2^ WA and f2g 2 ^. We 
compare the dynamics obtained by a numerical solution 
of the Floquet-Bloch-Rcdficld master equation using 
the exact eigcnstates of the Floquet Hamiltonian, the an- 
alytical Van Vleck-MRWA approach, Eqs. (gSJ) and 0S]), 
and the RWA. While the Van Vleck and RWA solutions 
relax incoherently to a stationary state, the numerical 
solution exhibits two full oscillations with VL m . As in the 
nondissipative case, the exact oscillation frequency f2 m 
is nonzero. For stronger damping those slow oscillations 



disappear. Both the numerical and Van Vleck oscillations 
show fast driving-induced oscillations which survive also 
in the stationary state. While for short time scale, Fig. 
IToT b). both approaches agree quite well, one finds that 
in the long time limit the amplitude of the fast oscilla- 
tions predicted by the analytical solution is smaller than 
the exact numerical one. Compared to the RWA solu- 
tion, where the equilibrium value is reached after longer 
time and the fast oscillations are averaged out, Van Vleck 
perturbation theory is clearly an improvement. 
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FIG. 17: (Color online) Driving- induced tunneling oscillations for the nondissipative case (k = 0). The survival probability 
P±-tl(t) is shown for A/ A = 3.0, co/A = 2.0, and e/A = 5.9011 (exact 3-photon resonance). Three approaches are compared: a 
complete numerical solution of the Floquet Hamiltonian, the second-order Van Vleck approach, and the RWA approach. For the 
first two approaches, complete population inversion is predicted, and for the Van Vleck dynamics we find the main oscillation 
frequency O3 = A| J-^{A/ijj)\. Besides, the numerical and Van Vleck approaches exhibit small driving-induced oscillations, see 
especially figure (b). The RWA predicts a strongly shifted oscillation frequency. Moreover, population inversion is incomplete. 
For further comparison, the RWA and the Van Vleck approaches are shown for a slightly shifted external frequency, u)/ A = 1.9. 
The dynamics in this case are almost completely localized. 




FIG. 18: (Color online) Driving-induced tunneling oscillations for the dissipative case (k = 0.01 and h/3A = 10). Remaining 
parameters are the same as in Fig. 1171 The numerical solution of the master equation ([38]), the analytical Van Vleck-MRWA 
solution, and the RWA solution for oj/A = 2.0 are compared. Good agreement between the numerical and the Van Vleck 
solutions can be observed on both long (a) and short (b) time scales. Within the RWA, not only are the main frequency and 
amplitude changed and the driving-induced oscillations missed, but also the equilibrium value lies far above the Van Vleck 
prediction. The RWA and the Van Vleck-MRWA solution for u/ A = 1.9 show an almost incoherent decay, and their long time 
limits differ strongly from the corresponding ones for uj/A — 2.0. 



D. Driving-induced tunneling oscillations 

An effect contrary to the CDT are driving-induced tun- 
neling oscillations (DITO). It has been predicted in [fUl] 
and experimentally shown in Q that for a high static en- 
ergy bias, e 3> A, and for high driving frequency, cj ^> A, 
coherent oscillations with frequency A\J- m (A/uj)\ and 
large amplitude are induced if e w mui. The DITO are 
often also named Rabi oscillations even though in the 



original problem of Rabi [78| a circularly polarized driv- 
ing field couples to the TLS. As a consequence, the ob- 
tained frequency of the oscillations is linear in A. 
In this section, we are going to investigate the effect in the 
regime of moderate energy bias and driving frequency. 
First, we examine again the nondissipative case (n = 0), 
see Fig. [T7] As parameters, we choose a moderate driv- 
ing amplitude and frequency: A/ A = 3.0 and to/ A = 2.0. 
For an exact 3-photon resonance, condition (|21[) must be 
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fulfilled and thus e/A « 5.9011. Notice that the RWA 
resonance condition (J7J, £ — 6cj, is only valid in the case 
of high frequencies u> ^> A. With condition (|21|) used, the 
Van Vleck approach results in the oscillation frequency 
fl { 3 2) = A\J- 3 (A/u>)\, and for times t = (2n + l)7r/fi 3 2) 
with n = 0,1,2,3..., one finds complete population in- 
version, see the Van Vleck graph in Fig. UTT a). Further- 
more, in Fig. 117( b). one can nicely see the modifications 
resulting from the external driving: three small oscilla- 
tions corresponding to a 3-photon resonance. The exact 
numerical solution shows a slightly shifted main oscilla- 
tion frequency The RWA approach exhibits the os- 
cillation frequency ftf WA , which is strongly out of phase 
compared to the numerical and Van Vleck one, and also 
has a smaller amplitude, so that a complete population 
inversion is not reached. When changing the driving fre- 
quency to slightly out of resonance, the driving-induced 
tunneling oscillations are strongly suppressed, and the 
system is almost completely localized in the initial state. 

This behavior originates - contrary to the CDT - not in 

(2) 

a zero oscillation frequency f2 3 but rather in a vanish- 
(2) 

ing amplitude of the fi 3 oscillation. Also the RWA at 
u = 1.9 A is suppressed. 

In Fig. 1181 we consider the influence of the environment. 
At exact resonance, we observe within the numerical solu- 
tion and the Van Vleck-MRWA approach coherent oscil- 
lations decaying to a stationary equilibrium value. Before 
reaching the equilibrium value, the dynamics are domi- 
nated by the frequency ^3, while in the long time limit 
the coherent oscillations die out; faster ones with the 
driving frequency u> and its higher harmonics around a 
static equilibrium value are found. The agreement be- 
tween the numerical and analytical calculations is quite 
good. Also, in the RWA approach, the coherent oscil- 
lation of frequency f2g ,WA dies out to a stationary state. 
However, apart from the frequency shift already observed 
in the nondissipative case and the smaller amplitude, the 
equilibrium value also differs strongly from the one ob- 
tained within the Van Vleck solution. Furthermore, since 
fast oscillations are completely neglected, the stationary 
state is constant. Considering the Van Vleck solution for 
a slightly shifted driving frequency, uj/A = 1.9, we no- 
tice an almost incoherent decay to an equilibrium value 
which is much lower than the one of the dynamics with 
uj/A = 2.0. Thus, dissipation leads here to an almost 
complete inversion of the population. 
We observe that our analytical methods are also able 
to recover the findings for the population difference in 
Chapter 3.2 of [|| in the high-frequency limit (u ^S> A) 
and even can reproduce the small modulations which are 
found there by a numerical treatment of the dynamics. 
Furthermore, we are able to go beyond the assumption 
of a high driving frequency. 



IV. CONCLUSIONS 

In conclusion, we discussed the dynamics of the spin- 
boson system exposed to an external ac driving. Assum- 
ing weak coupling between TLS and bath, we arrived at a 
closed analytical expression for the time evolution of the 
system. Our results are at the same time valid for the 
whole range of the driving amplitude A and for moderate 
to high driving frequencies ui, see discussion in Sec. Ill Dl 
In contrast to the NIBA, we are able to treat both an un- 
biased and a biased TLS for low temperatures and weak 
damping. Indeed, besides the Born-Markov approxima- 
tion, the only further simplifications we used solving the 
time-dependent Hamiltonian are the moderate rotating- 
wave approximation in Sec. IIIIBI and the expansion in 
the dressed tunneling matrix element A„ = J n (A/u>)A, 
with A the bare tunneling coupling and J n the nth-order 
Bessel function, using Van Vleck perturbation theory. In 
the vicinity of an m-photon resonance, the latter is also 
justified, as shown in Sec. Ill Dt for moderate driving fre- 
quencies as long as condition (f2"5"j) is valid. We found cor- 
rections to the renormalized Rabi frequency fi^ WA [Eq- 
(fTT|) ] also leading to a shift of the resonance condition 
for an m-photon resonance [Eq. (|2 1 [) ] - The so-calculated 
quasienergy spectrum is in very good agreement with re- 
sults found by a numerical diagonalization of the Floquct 
Hamiltonian for all values of the static bias e. Upon in- 
vestigation of the survival probability Pj-^t), we could 
recover the shifted oscillation frequency reported already 
in [59|. We included also the second-order modifications 
to the Floquet states in our calculation, which account 
for the higher harmonics induced through the external 
driving and lead to fast oscillations in P|_^(t), see Figs. 
G2EE3 andDJ 

By adding a thermal bath to the TLS, we examined in 
Sec. IIIII the dissipative dynamics of the system. In Sec. 
IIII Al we visualized the good agreement between our an- 
alytical formulas for the position matrix elements and 
a numerical calculation even for low driving amplitudes. 
This turned out to be essential to arriving at a phys- 
ically realistic result for the relaxation and dephasing 
rates given in Sec. IIII Bl Comparing RWA to Van Vleck 
results, we found strong deviations and even unphysi- 
cal predictions for the former one at low driving ampli- 
tudes. We remark that our rates agree very well with 
the zero-temperature results derived recently in [3^] via 
the dressed state approach. In this work, a charge qubit 
is strongly driven by a microwave field and connected to 
a dc SQUID. The higher order corrections to the rates 
prove to be essential to correctly reflecting the physical 
findings in this experiment. From this we are encouraged 
that our results provide a realistic picture of relaxation 
and dephasing processes in a driven two-level system and, 
due to the generality of the model, are of interest to a 
wide range of physical applications. 
In Sec. IIII CI we performed a detailed analysis of the TLS 
at an exact 3-photon resonance and for a vanishing third- 
order Bessel function, which is known to lead to coher- 
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cnt destruction of tunneling in the high- frequency limit. 
For moderate driving frequencies, we found second-order 
modifications to the RWA solution. While the latter pre- 
dicts a complete localization of the TLS in the initial 
state, the Van Vleck solution shows that driving-induced 
oscillations survive. Furthermore, for the dissipative case 
we found an incoherent decay to a quasistationary value. 
In Sec. MI D[ we examined an effect opposite to the 
coherent destruction of tunneling: for an appropriately 
chosen driving amplitude, coherent tunneling oscillations 
with frequency TO , Eq. ([20]) . can be observed at an 
m-photon resonance. By slightly changing the driving 
frequency out of resonance, these oscillations are almost 
completely suppressed and the system shows an incoher- 
ent behavior. 
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is defined as 

«•!■» = k[ dt{ ' [) - (A5) 

Considering a T-periodic state vector \u a ^ n {t)) living in 
the spatial Hilbcrt space %, we can write it in a Fourier 
series and thus expand it in basis functions of 7~: 

K»(*)> = e-^'Mt)) = Y, e- Uut \u£- l) ), (A6) 

i 

where \u^} are the time independent Fourier coeffi- 
cients. In the composite Hilbert space H <E> T, we define 
the state 

|u„, n >>=$>a"°>®IO- (A7) 
I 

Through the expansion of the Hilbert space, we can now 
treat the time-dependent problem (|A1|) like a time inde- 
pendent one. 



Appendix A: Floquet theory 

The Floquet theorem states that the Schrodinger equa- 
tion 

mm)) = H(t)\m) (ai) 

for a Hamiltonian being periodic in time is solved by 

\Mt)) = \u a (t))e-^l\ (A2) 

where \u a (t)) = \u a (t + T)) and we assume that the oscil- 
lation period of H(t) is T = 2ir/u>. The quasienergies e a 
are obtained as eigenvalues of the Floquet Hamiltonian 

H\u a (t)) = [H(t) - ihd t ]\u a (t)) = e a |u a (i)). (A3) 

Note that \u a . n (t)) = exp(— muj£)\u a (t)) yields a solu- 
tion of (|A1|) physically identical to (|A2[) but with the 
shifted quasienergy e an = e a — hnu. Furthermore, 
£a = £a,o and \u a (t)) = |u a ,o(i))- Thus, it will be suffi- 
cient just to examine the set of eigenvalues {s a ,n} with 
— haj/2 < S a ,n < hui/2. 

We introduce the Hilbert space T of the T-periodic func- 
tions, with the inner product defined as 

{f,9) = ^f Q dtf*(t)g(t). (A4) 

The functions Lp n (t) = exp(— inuit) build an orthonor- 
mal and complete basis set of T [79| , where we further 
define for a basis-independent notation the state vectors 
\n) with ip n (t) = (t\n). The scalar product in the ex- 
tended Hilbcrt space % ® T of the Floquet Hamiltonian 



Appendix B: Van Vleck perturbation theory 

Here, we give the matrix elements of the effec- 
tive Hamiltonian Ji c s and the transformation matrix 
exp(±iS) = 1 ± iS^I ± iS ( V + iiSWitfC 1 ) to second- 
order in A expressed in the eigenstatcs ([3]) of the unper- 
turbed Hamiltonian. For first order in A, the Hamilto- 
nian has the same shape as within the RWA. Its elements 
are @, ^ 

««tA,nl W 2K/W» = ^nKu (Bl) 
««t,n|w2l<i» = -|A„_ ; 5 n _ iim . (B2) 

The elements of the transformation matrix are 

«4JiS (1) K.*» = \ + A T l n (l - S '-n,m), (B3) 

1 2 e + u)(n — I) 

*' 1 L £ + u!{n — I) 

The Kronecker 5 comes from the fact that iS 1 ^ vanishes 
between almost degenerate states by construction. For 
the second-order elements, we find 

««?^ b i«2k /w » = t5 E jrL (bb) 

l^ — vri 

(Knl^ffKn+J) = ««!,n+ml^K.n» = 0- (B6) 

The expression for the transformation matrix already be- 
comes more evolved: 
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(<4Jis< 2) K,,->} 



i 



4(n — j)uj 



E 

k^j-\-m 



^■n—k^j — k 



i 



i 



^n—j—m^—m 



^j—n—m^—m 



+ (ft — j — m)uj e + (j — n — m)uj 



e + (n — fc)a; e + (j — fc)o; 

(i - <w, 



(B7) 



(«„|tf (a) i«!j» 



1 f \ -> Afc_„Afc_j- 

kjtj-m 



1 



Aj—m—n A_ m 



A n — TTj — j A_ m 



-e + (n — j + m)w — e + (j — n + m)w 



e + (n — fc)w — e + (j — k)uj 



(B8) 



(«n|i5 (2) |«w» = ((«wl^ p) K,»)) = °- 



(B9) 



By applying the transformation now on the eigenstates of the effective Hamiltonian |<&?? n )), see Sec. Ill CI we get the 
eigenstates of the Floquet Hamiltonian "Htls to first order in A, 



2 e + jco 



sign (A_ m ) cos — \ul j+n+m )) - sin — lu^^-)) > , (BIO) 



l<0> =!*£»» + \ E TTfc {^(A- m )sin^|4 j+n )) + cos ^ !<„_„,_,-))} • (Bll) 



And to second-order, 



e 



e, 



l*L 2) „)) - |*™» + £ sin -fK i+n » + sign (A_ m ) cos -f |<_ i+n+ro )) 



x < 



A; 



A-j-m 



k^—m j^ — m 



e + (j — m)cj c — (j + m)u; 

AfcA,- r . e 



E 



p^-j—m 
p^—m 



1 Aj_|_pAp 
2 



1 



1 



£ + (j + P)u e + puj_ 

e 



sin^K.Jfe+n-a)) + si g n ( A -m)cOS^|4 1 fe+n+m-i)> f 



(B12) 



l*S»» = |^)>-^{cos^|< 3 . +n _J)-sign(A_ m )sin^|<_ j+n >)} 




Aj— A_ m 



A_j_ TO A_ m 



e + (j — m)cj e — (j + m)w 
A fc Aj f e 



E 



1 Aj+pAp 



1 



1 



(e + ku)(s + juj) 



4juj 2 

p^—rn 

-y 1 K.fc+n-j-m)) ~ si g' n (A- m ) sin ^ i u 



£ + (j + p) w e + puj 

e 



(B13) 



Appendix C: Calculation of the dynamics p(t) of the TLS, fulfilling the condition that p(0) = | t) I - 

By diagonalization of the Floquet matrix Q or by solving 

To calculate the survival probability of the system, 
P 7 _v 7 (£), where 7 =1~, J,, we start with the density matrix 
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the master equation (|38[) . we obtain the density matrix 
in energy basis with the matrix elements 



so that 



p a p(t) = ($ a (t)\p(t)\*(,(t)) a,P = ±. 



(CI) 



Using that p (t) + p++{t) = 1 and p~+(t) = p* + _{t), 

we get 

P^{t) = 2Re{( 7 |$_(t))($ + (t)| 7 >P-+(t)} 

+ |< 7 |<M*)>| 2 + (|<7lM*))| a - l(7l*+W>| 2 ) P-(t). 

(C2) 

The corresponding transition probability is just 
P 7 ->i(i) = 1 - P 7 _> 7 (t), where S ^ -f. From jC2]), we 
see that we have to calculate (7|$ Q (£)). We use the 
periodicity in time and express it in a Fourier series: 



( 7 |* a (*))=X;<7l*? ) )«p(iftwt) > 



(C3) 



with 



(7l^ fc) ) = ~ / rfte- ifc " 4 (7|<f> Q (i)) 



TJo 



dt(>y\* a , k (t)) = ((7|$a,fc», (C4) 



where we used in the last step the definition for the in- 
ner product of the extended Hilbcrt space, Eq. (|A5|) . 
and defined | 7 )) = | 7 )|0). By this we establish a connec- 
tion between the Floquct states and the time-dependent 
Hilbcrt state: 



{ 1 \<£ a (t))=J2((l\<S> a , k )W kuJt - 



(C5) 



-(*) 



-(0), 



and 



p_+(t) = p_+(0) exp [i (mcj + t 



(C7) 



(C8) 



where we used the general expression for the quasiencr- 
gies at an m-photon resonance found in Sec. HH The 
starting conditions are calculated through 



Pafl(P) = <$ a (0)|7><7l**(0)>. 
Combing this, one gets 



(C9) 



P^(t) = (|( 7 |$_ (t))| 2 - |< 7 |<M*)>| 2 ) l(7l ( f>-(0)>| 2 
+ |(7|<f> + W)| 2 + 2Re|( 7 |<i>_(t)>(1> + (t)l7> 
x ($_(0)|7)(7|$ + (0))c i ( m "+°-) t l. (CIO) 



a. flWOl survival probability 

Using in this general expression the eigenstates (fT2")l 
and (|13p . we arrive at the survival probability in the 
RWA, 

(Gil) 



1. Survival probability in the nondissipative case 



For the Hamiltonian of the nondissipative TLS, Eq. 
([1]), the master equation is simply 



Pap{t) = -K £ a _ £p)Pa.p(t), 



(C6) 



6. Van Vleck survival probability 

To get the survival probability to second-order in A, 
we use (|B~T0|) - (|B"T3)) in flC5j} and obtain 



(I \^\t)) = exp (-iA- s inu>t) |-sign(A_ m )cos^e- iro -* - I sin ^(f) 

+ sign (A_ m ) cos £p e - imwt [B(t) + C(i)] J , (C12) 



(I |$i 2) (i)> = exp (-i^ mnuitj |-sign (A_ ra ) sin ^ + 1 cos ^ e imw U(i) 

+ sign (A_ m ) sin % [B(t) - C*(t)]\, (C13) 
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where we defined 



Anuit 



n^—m 



e + row ' 



(C14) 



B(t) = Ve inu ' — 
4nw 



A n _ m A_ 



- m — n — m 



s + (n — m)u e — (n + m)w 



(C15) 



c(t) = y y ApAp+ " e ! "" f 

w L^, 8nw 

n^O p^ — m 
p^—n—m 



s + pco e + (p + n.)u; 



(C16) 



Using those expressions in (|C10j) . we obtain the survival probability 

P^(t) = P™f W + P^(t) + P<\(t). 



(C17) 



We distinguish three different parts. The first one corresponds to the averaged second-order Van Vlcck approach: 

(C18) 



p™f (t) = cos 2 (n£>|) +cos 2 e m sin 2 (n£> 



Then we have additional contributions from Van Vlcck perturbation theory to first order in A: 

A. 



Pl%{t) = -isign(A m )sine m sinfi^i ^ 



n^ — m 



e + nui 



sin[(ro + m)ujt\ 



+ sign (A m ) sin 26 m sin 2 ( - 



n^—m 



(C19) 



And finally the second-order part: 



P^it) = -sign (A m ) sin 26 m sin 2 ( f2^| ) J-U(O) 2 - 2[B(0) + C(0) cos 6 



sin 2 8 m sin 2 fll 2) 



3 E 



4 \ e + row 



sin(n + m)wi J + — I 



£ + row 



cos(ro + m)ujt 



+ I 1,4(0) J \A{QY 2[B(0) + C(0) cos 9, 



+ 1,4(0) - J 1^(0) 2 - 2[S(0) + C(0) cos 0, 



E 



n^—m 



£ + nw 



cos(ro + m)wt cos 2 (f2^ ) 2 - ) -J + cos0 m sin(ro + m)ujt smilfft — cos(ro + m)u)t cos 2 m sin 2 (^m^ 



^)l)-cos 2 e m sin 2 te| 



1 cos 2 te)l)+cos 2 e m sin 2 k( 2 )l 



E- 



2row 



A 7i — m A_ m A_ m _ n A_ 



E E 



e + (ro — rro)w e — (ro + m)w 
AjAfe 



cos nwi 



(e + fcw)(e + jw) 



cos(j — k)ut — f cos 6, ; 



E E 

p^ — n — rn 



+ £ + (p + n)c(j 



j^L—m k^ — m 

cos B m cos nujt — 2d 



^)l)+cose m sin 2 (n£)i 



(C20) 
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Here, for e 7^ muj, we have 



ia 2 - 26- 2c cos 9^ 



4ra - ^ba 2 + 26 + 2c 2 - 2c(i a 2 - 36) cos 9 2 + 2c 2 cos 29 2 
16 4 2 



— -a l -a 2 — 26 — c — 3c cos 9 2 ) \ -a 2 — 26 — c — 3c cos 9 2 
2 V 4 Tn / V 4 



and 



/ = T 



a 2 - 26- 2c cos 91 



ca 2 - 26c - 2c 2 cos Qt 



with a = A(0)/A, 6 = £?(0)/A 2 , and c = C(0)/A 2 ; while in the case e = muj, the definitions 

d = B(0), f = -B(0) 

have to be used. 
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